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TECHNICAL NOTE 2129 

METHOD OF CALCULATING THE LATERAL MOTIONS OF AIRCRAFT 
BASED ON THE LAPLACE TRANSFORM 



By Harry E . Murray and Frederick C • Grant 


SUMMARY 


The lateral motions of aircraft are obtained by means of the 
Laplace transform which gives solutions expressed in terms of elementary 
functions for the free and forced motions. These equations permit the 
calculation of the free motion of an aircraft following any initial 
condition or the forced motion following the application of constant 
external forces and moments. These forced motions can be used to obtain 
by means of Duhamel's integral the response to any arbitrary forcing 
function. All the classical stability concepts can be deduced from 
these same solution equations largely by inspection. These equations 
for the lateral motion are applied to the calculation of the lateral 
stability of a specific airplane and to the calculation of certain of 
its free and forced motions. 


INTRODUCTION 


The lateral motions of aircraft are represented by three simultaneous 
differential equations which are generally assumed to be linear. The 
fundamental problem of lateral dynamics involves the solution of these 
differential equations in terms of the aerodynamic and mass parameters 
of the airplane. The solutions can then be used to obtain numerically 
the motion of the airplane as a function of time. 

The recent application of the Laplace transform to the solution of 
systems of linear differential equations permits a more general analysis 
of the problem of airplane motion than that of reference 1, which is 
based upon Heaviside's operational calculus. Heaviside's operational 
calculus permits a calculation of the forced motion, which is the motion 
following the application of external forces and moments. The Laplace 
transform permits these same calculations and also permits the direct 
calculation of the free motion, which is the motion following finite 
initial values of the variables and their first derivatives in the 




2 


NACA TN 2129 


absence of externally applied forces and moments. This calculation 
cannot be made by use of Heaviside's operational calculus. The Laplace 
transform solutions, which include both the free and forced motions, may 
be written in a closed form from which a ll the classical stability con- 
cepts can be deduced largely by inspection. The form of the equations 
of motions of the airplane is independent of such aerodynamic parameters 
as Reynolds number and Mach number, and these parameters enter the equa- 
tions only as they effect the values of the aerodynamic constants or 
stability derivatives appearing in the equations. The values of the 
stability derivatives must be obtained by actual measurements during 
physical tests or from aerodynamic theory before motion calculations can 
be attempted. 

Investigations of some of the possibilities of applying the Laplace 
transform to the study of aircraft motion have been reported in refer- 
ences 2 and 3* and in two British reports, one by K. Mitchell, the other 
by J. Watham and E. Priestley. The British papers do not give final 
equations in a form sui table for calculation purposes. The analysis of 
reference 2 closely parallels that of the present paper until the point 
of taking the inverse Laplace transform is reached. At this point, 
reference 2 indicates that the inverse Laplace transform can be taken 
either by means of the relatively simple partial-fraction expansion 
(used in the present paper) or the more complicated inversion theorem 
of the Laplace transform. Neither approach in reference 2 is carried 
to the point of final equations containing only elementary functions and 
in a form particularly suited for computation. A solution similar to 
that of the present paper is indicated in reference 3. Only the form 
of the analysis is shown in. reference 3 > however, and all the details 
necessary for practical applications have not been carried out. 

The present paper presents an analysis based on the representation 
of the lateral, motion of an aircraft by differential equations. The 
results of the analysis are solutions in closed form expressing the free 
and forced motions in terms of elementary functions. These equations 
permit the calculation of the free motion of an aircraft following any 
initial condition or the forced motion following the application of 
constant external forces and moments. These forced motions can be used 
to obtain, by means of Duhamel's integral, the response to any arbitrary 
forcing function as shown in references 4 and 5- The solutions are 
readily adaptable to calculation by digital-type calculating machines 
and the calculation is an arithmetical process requiring no knowledge of 
the theory of the Laplace transform. The solution equations of motion 
have been applied on an automatic calculating machine to the calculation 
of the lateral stability of a specific airplane and to the calculation 
of certain free and forced motions as illustrative examples. 



COEFFICIENTS AND SYMBOLS 


trim lift coefficient (W cos 7/q£3) 

rolling -moment coefficient (L/qSb) 

yawing-moment coefficient (N/qSb) 

lateral-force coefficient (l/qS) 

airplane weight, pounds 

rolling moment 

pitching moment 

yawing moment 

lateral force 

aileron hinge moment 

elevator hinge moment 

rudder hinge moment 

dynamic pressure (pV2/2) 

wing area, square feet 

wing span, feet 

inclination of flight path to. horizontal (positive 
in climb), degrees 

angle of attack, degrees 

angle of pitch, degrees 

mass density of air, slugs per cubic foot 
free-stream velocity, feet per second 
airplane mass, slugs (w/g) 

acceleration due to gravity, feet per second per 
second 
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8b 


nondimensional time (tv/b) 


t 



4 

Pb 


0 

* 


t im e, seconds 


inclination of principal longitudinal axis of inertia 
(positive for axis above flight path at nose), 
degrees 

airplane relative-density factor (m/pSb) 


angle of bank, radians 



angle of yaw or azimuth, radians 



r dt) 


P 


rolling velocity about stability X-axis, radianB per 
second 


r 


yawing velocity about stability Z-axis, radians per 
second 


P 


angle of sideslip, radians 



rolling-moment coefficient of forcing- function couple 
in roll 


^n 


c 


C Y c 

P,P' 

T l/2^ T l/2* 


N l/2> N l/2' 

^a 

5r 


yawing -moment coefficient of forcing-function couple 
in yaw 

lateral— force coefficient of lateral forcing function 

periods of oscillatory modes, seconds 

times to damp to half-amplitude of oscillatory modes, 
seconds 

cycles to damp to half -amplitude of oscillatory modes 

aileron deflection, degrees 
rudder deflection, degrees 
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&e 

K X 


K Z 


K xz 


kx 0 

% 


a 

A 

X 2 , Xg > 


R a } 


V 


dq^i 

da 


R 1 


I 


1 


R 1 


t 


elevator deflection, degrees 

nondimens ional radius of gyration about stability 
X-axis 




nondimensional radius of gyration about stability 


Z-axis 


7 q 2 cos 2 t] + Kx^sinO 


nondimensional product of inertia between stability 
X- and Z-axes ^KZq 2 - Kxq^ sin T] cob 

non dim ensional radius of gyration about principal 
X-axis (kXo/b) 

nondimensional radius of gyration about principal 
Z-axxs (kz 0 f°) 

radius of gyration about principal X-axiB, feet 
radius of gyration about principal Z-axis, feet 
Laplace transform of s^ 
stability quart! c 
roots of A = 0 


polynomials in a 


roots of q^ = 0 
real part of and 


when 

X 3 

and 


Eire 

complex conjugates 
imaginary pant of X .3 

and 


when 

x 3 

and 


are complex conjugates 
real part of and ^ 

wben 

X 1 

and 

*2 

are 

complex conjugates 
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I]_' Imaginary part of X-]_ and Xg when X^ and Xg 

are complex conjugates 

A,B,C,D,E coefficients of stability quartic 

R Routh' s discriminant 


A 1> A 2 ,A 3 ,A V A 5 ,A 6 

amplitude coefficients for 

0 





Bi,b 2 ,B3,B4,B5,B6 

amplitude coefficients for 

* 





01,02,03,01^,05 

amplitude coefficients for 

P 





% 

real part of A 3 and A 4 
complex conjugates 

when 

x 3 

and 

X4 

are 

*A 

imaginary part of A 3 and 
are complex conjugates 

A k 

when 

x 3 

and 

X4 

*b 

real part of B 3 and B 4 
complex conjugates 

« 

when 

x 3 

and 

X4 

are 

Ib 

imaginary part of Bg and 
are complex conjugates 

B 4 

when 

x 3 

and 

H 

Rc 

real part of C 3 and Clj. 
complex conjugates 

when 

X 3 

and 

Xlj. 

are 

*C 

Imaginary part of Cg and 
are complex conjugates 

C 4 

when 

x 3 

and 

X 4 

V 

real part of A]_ and Ag 
complex conjugates 

when 

X 1 

and 

Ag 

are 

V 

imaginary part of A^_ and 
are complex conjugates 

A 2 

when 

X 1 

and 

Xg 

V 

real part of B^_ and Bg 
complex conjugates 

when 

X 1 

and 

X2 

are 

Ib’ 

imaginary part of B^_ and 

Bg' 

when 

X 1 

and 

Xg 


are complex conjugates 



real part of and 

complex conjugates 


C 2 when and 


^2 are 


Imaginary part of and C 2 when X.-j_ and X2 

are complex conjugates 

amplitude coefficient for 0 oscillation corre- 
sponding to complex conjugate roots X3 and X4. 

(2 ^R A a + Ia 2 ) 

amplitude coefficient for \jr oscillation corre- 
sponding to complex conjugate roots 7.^ and X4 

(2 \/r b 2 + I B 2) 

amplitude coefficient for p oscillation corre- 
sponding to complex conjugate roots X3 and X4 

(2 yfe c 2 + ic 2 ) 

amplitude coefficient for 0 oscillation corre- 
sponding to complex conjugate roots X^ and X2 

(2 yV 2 + i A ' 2 ) 

amplitude coefficient for t oscillation corre- 
sponding to complex conjugate roots X x and X2 

(2\Ab’ 2 + I B ’ 2 ) 

amplitude coefficient for p oscillation comer- 
sponding to complex conjugate roots X^ and X2 

( 2 \Ac' 2 + i c ' 2 ) 

phase angle for 0 oscillation comesponding to 
conjugate complex roots X3 and X4, radians 



phase a n gle for oscillation comesponding to 
conjugate complex roots X3 and X4, radians 

/I -1 ib\ 
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“A 


I 


(Ob’ 


°c' 



phase angle for 3 oscillation corresponding to 
conjugate complex roots X3 and X4, radians 



phase angle for 0 oscillation corresponding to 
conjugate complex roots and X^d radians 

^fcan - -'- 

phase angle for oscillation corresponding to 
conjugate complex roots and Xq, radians 

^tan - l 

phase angle for 0 oscillation corresponding to 
conjugate complex roots X3. and X2, radians 
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i _ 

' n p ■ ar 



coefficients appearing in numerator terms of ampli- 
tude coefficients for 0 

coefficients appearing in numerator terms of ampli- 
tude coefficients for ijr 

coefficients appearing in numerator terms of ampli- 
tude coefficients for p 


initial value 
transformed variable 


dCv 

% = ST 


a o* a i> a 2J a 3J a ip a 5 




C 0d» c l-> c 2> c 3> c 4 


Sub scripts: 

0 

cr 


ANALYSIS 


The linear equations of motion, referred to the axis system shown 
in figure 1 and representing the lateral motion of an airplane are 



10 


NACA TN 2129 


2 p b K x 2 D b 2 0 - + SMbKxzPb 2 * " TPlyPb* ~ C lJ ~ C l r _ = 0 


2 4 r 


2^X2^^. - ^np^ + 21^^%% - ^Dfct - C^P - = 

-io^Dbtf - + 2^+ _ c L tan 7 * - ic^t - C Yr P + 


/CD 


2 p b D b p - C Y = 0 


The terms Cj , C^, and Cy c are forcing functions which represent 

disturbances imposed upon the state of motion of the airplane hy control 
movement or atmospheric turbulence. These terms , in general, are arbi- 
trary functions of time, but for the purpose of this analysis, they are 
considered to be constants applied at zero time. After a solution has 
been obtained in terms of constant forcing quantities this solution can 
be used to obtain a new solution for an arbitrary forcing function by 
Duhamel's integral as explained in references 4 and 5. 


Transformation of Equations 

When the Laplace transform is applied (reference 6, p. 8), the 
transformed equations become after multiplying through by 0 


(2hbK x 2 03 _ |cz p a2)0 CT + ( 2 ^x 203 - + (-C^^ = r ± 


r l = (2hbK X 2 0 2 - £Clpa)0 O + (2^X202 - ^0)^ + 
( 2 Mb K X 2cT ] (^) 0 + { 2 ^JZ a ) (^0 + C l c 


V ( 2 a) 
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( 2 ^xZ ct3 " | C V 2 )^ 0 + ( 2 ^Z 2a3 " 2°“/)^ + (" C V) Pa = r 2 
r 2 = (2^bKxZ^ 2 - ^n p ^ 0 + ^h^a 2 - |c nr a) * 0 + > (2b) 

(2PbKxZo) (M) 0 + (2MbKz 2 ^)(l>bV)o + C n c ^ 


(_ic Y ^ CT 2 _ c L a)0 a + (2^a 2 - ic^a 2 - C L tan 7a) ^ + 
( 2 W 2 - CY p a)p 0 = r 3 

r 3 = (-F C Y p a )^0 + ( 2ti b CT " ^Yr 0 )^ + ( 2ti b°)P0 + C x c 







Solution of Transformed Equations 

After equations (2) are solved by determinants, the expression for 

fa ls 

^ + a^a^" + &2.a^ + a^a 2 + aj^a + a^ 

<f^A 


where 


A = Aa^ + Bff3 + ca 2 + Da + E 


(*) 
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and the constants axe given "by 


ao = 0 o l 8 ^ 3 


K x 2 K JZ 

%Z H* 


/ 


a l 


= 00 


2 ^ 


\ 


Kx 2 K XZ c l 


K- 


XZ A Z 


x 
L x 
-2Cy 


Cn T 


\ 


- 2\*b 


C Zp K XZ 

% Kz 2 


(Db0) O 8Mb 3 


K x 2 K XZ 
KXZ «Z 2 


ag - 0 O 



% 

% 

K XZ 


% 

Kx 2 

°V 


Cz p Kj2 
C np K XZ 


Cnp 

Cn P 

*Z 2 

+ Pb 

c np 

K xz 

Cn r 

- >m b 2 

\ 

% 

Cy p 

0 


ccl 

& 

0 

C *x 



\ 




% % 

C n-n C n r 


(Db 0 )o P* 


7 


(l>bt)o^ M b 2 


c Z r K XZ 



Cz c %Z 
Civ KZ 2 
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1% CZ 3 


-2 tan 7C l 


~ I C np C n p °n r + 

% % / 


/ % ^ C * r 
(^) 0 k c np %Z C n r 

\ % 0 C Yr 


o C *B K X® 
- %tb 2 P . 

c np K XZ 


x / C 2p K XZ C i r 

/ Cj K XZ \ 

i 0 /2^ tan 7 Cl P 2 I ^ Cn p K Z 2 Cn r “ 

\ S Kz / <*„ o C r 


C ^p K xz 

c np K Z 2 


c z r % 

Cn r C nc 


/ \ / C Zp K xz Cz c 

+ Pof-I^b ^ lr ] + 2Mb Cbp K Z 2 + 

\ ** °*rj \ % 0 <*„ 



Ik 
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a 4 



1 

2 


tan 7C l 




(Db^) 0 2Mb 


tan 7 Cl 



K X 2 

K xz 


+ 


*0 


(- 


1 

2 


tan 7 Cl 


% 

c np 



+ CMOop-'b tan 7 C l 


l G Z c 


yn P 




C Z C 2 
t c 

Cn c C n r 
C Y C C Y r 



*5 


tan 7 Cl 


i% 

c np 


C ^c 

Cn c 


■A 



| K XZ 


K xz 

K z 2 



K x 2 

K xz 

C Z r 



OJ 

df 

OJ 

II 

pq 

K xz 

K z 2 

Cn r 

OJ 

£ 

CVJ 

1 

c z p K xz 


0 

1 . 



S K z 2 



Cj P % C * r 


0 K x 2 K xz 

o 

II 

cf 

0 Kje K z a 

+ »*h 

c np c np c nr 


^p % °*r 


% C Yp c Y r 


+ §**b 


% C *i 

Pn-n c n T 


- W 


Cz p K X C 


c n p K xz 
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D = 2Mt,CL 


E = - 


l r 

*2 C L 


% 

Cnp 

0 


"np 

0 


Oi 

** 

K xz 

K XZ 

OJ 

w 3 

1 

tan ' 

C7 

6 P 

C2 r 

Cn 

n P 

c n r 

1 

tan 7 


c h c h 

Cn p % 
% 


Cn- 


C Yt 


+ Vb 


C 7 C 7, 

c np c n p 


The expreBslon for y a is 


•ir = 
v cr 


hpcr^ + + bgcr^ + bgg 2 + blj.g + b^ 


g 2 A 


( 5 ) 


where the constants are given hy 


b p = *6 Rb 3 




K. 


XZ 




1 

K X 2 

K xz 

C *r 

h x = + 0 

2Pb 2 

K xz 

Y 2 

K z 

Cn r 


\ 

0 

1 

" 2Cy P 



0 >bt) o 8pb : 


K X 2 

%Z 
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t»2 = 


/ 

% ^ 

+ 

-<v 

o 

/ 

[-2Mb 2 


Mb 

l 

Cn p K XZ 

iy ' 

V 


0 K x 2 

C np Cnp 

;°Yp C Yp C ^r 


Mb 


Clp C Zp C lr 

o Kxz K z 2 

C Y Cy C Y 
P P r 


W 


ca 

O 

+ ^Mb 

% C *r 


+ (Db^)o 

( 

-2Mb 2 

Kx 2 Cz p Kxz 

Kxz C]^ K z 2 

|Cnp K XZ 

e± 

c n p C^. 

1 


l 

1 -2Cip 0 



| c l c Kx 2 



*bo — 


00 ^" 2 ^ c L 


% K x" 

C np K XZ 


(Bb0) o 


Mb 


\ 


Clp C Ip Kx 2 




c n T) K XZ 


°Y p % 


0 



% C Zp Kxz 
Cnp C np K z 2 

1 

" !<- 

C 7 C 7 c 7 

ip ip 4 r 

Cnp C Qp Cn r 


+ (Dbt) 0 

f 

Mb 

Clp C Zp Kxz \ 
Cnp C np K z j 


0 2C L 1 


c Yp C Yp CY r 

J 


V 

% % ° j 


Po Mb 


c 7 c z 

ip 6 P 
c np c n p 


/ 


2 Mb 


\ 


C 7 Ci ^ 


Cnp C nc Kjz, 

°Yp <lC Yc ° 


+ Mb 


C 7 C l-n 

t C p 
C nc CrL p 


\ 


/ 
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^4 = 0oi 





+ C^b^o 


l-2(Jt|CL 


% 

Cnp 


Kx 2 

K xz 


(Db*) 0 





+ 


|% 

lC np 

% 


c h °*c 

C Hp C n c 

% 




The expression for p a is 


CnCT^ + c 


Prr = 


_ ^0 


,a3 + c o 0 2 + c o 0 + 


0A 


where the constants are given by 
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C 1 = 00 


+ (Db0) o 2(ib2c Y 



% 2 Kxz 

K xz K z 2 


(Db'Jr) 0 (2p b 2^Y r - 4^) 


K x 2 K xz 

*XZ K Z 2 


^^y ■ 


K X* Kxz 

K xz K z 2 


/ 

C Z r Kx 2 


c 2 p K xz 

^ + (Db0) O 


MbC L 


- MbC L 



V 

c n r K XZ 


% V 



2 V 


K x 2 % K XZ 
K XZ c n p Kz 2 
0 2C t 1 


\ 


Kx 2 % Cz T{ 


Up C n r 


Kxz C 

0 Cy Cy 


+ ^0 ( Vb tan ? C L 


1 


°i r % 2 

Civ Kxz 


Pb tan 7 Cj, 


C l v K XZ 


Kn 


/ 


(Db*) 0 


f*b 


V 


K xz 

% 

C *r 

K z 2 

° n p 

Cn r 

0. 

° Y P 

Cy 

x r 


2 V 


K- 


X 


K- 


•XZ 


j n r 


A XZ 

K z 2 


-Mb 


1 -2 tan 7 Cl 0 

0 Kx 2 Kxz 


/ 


+ P0 &Mb 


C Z C 7 
fc p 4 r 

c n p c n r 


Cjic C np C nr 

Cy Cy Cy 
x c i p • L r 


- Mb 


C z c c Z p Cz r 
• 0 Kxz Kz 2 


jCy Cy Cy 
± C p X T 


- w 


K X 2 C Z 


c 


K XZ C nc 
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c 3 = 0 o[^l 


% C *r 

c n p C nr 


+ O^b^o 



KX 2 % 

C *r 

MbC L 

K XZ c n p 

Cn r 


0 1 

tan 7 y 





\ 

/ 

K XZ 

Cz p 

c 2 r 



Cl v 

c 2 r 

\ 

I 





tan 7 Cl 

c n p 

c n r 

+ (Db^) 0 

MbC L 

K Z 2 

c n 

n p 

Cn r 





/ 


0 

i 

tan 7 



C l c C Zp C Z r 

Cric C n p Cn r 

CY c CY p Cy r 


- 2MbC L 


C 2 c Kx 2 K XZ 

C nc Kxz Kz 2 

0 1 tan 7 


+ Mb 


% C *c 
Cnp C nc 




c 4 = 5GL 


-C “p r 

Cn p C n r 
0 1 tan 7 


All the determinants given in thiB paper are expanded in the 
appendix. 


In order to obtain the actual variables 0, i|r, and (3 from the 
transformed variables an inverse LaPlace transformation must be applied 
to fa> ^q; and Pa* The expressions for 0 a , i a , and 0 0 are of the 
form p a /qc where p 0 and are polynomials, the degree of 

being higher than that of p 0 . Reference 6, page 45, indicates that the 
inverse transform of a function of this type is (in terms of the variables 
used herein) 


L 

W 


Pq( 4n) ^ ln s b 

n=l 9q ' (^-n) 


This equation assumes all the roots X n of qq = 0 to be distinct. All 
roots of qc = 0 are distinct for 0 0 ; however, for 0 0 . and tg, 

% = 0 has double zero roots. (See equations ( 3 ), ( 5 ), and ( 6 ).) The 
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terms in tbs equations for 0 and ijr resulting from the two zero roots 
of = 0 are given according to reference 6, page k-9, by 


32(o) + ^(oK 

da 


where 

Q = ?2a2 
<3a 

Ttip analysis takes three forms depending upon the character of the 
nonzero roots of q^ = 0 which are the same as the roots of A = 0. 
Four real roots, two real roots plus a pair of conjugate complex roots, 
or two pairs of conjugate complex roots may exist. 


Four Real Roots 

The inverse Laplace transform of $ a is 


■fi = + A^ 3 ” + A 3 e X 3 8 » + + A 5 b„ + Ag 


( 7 ) 


where the constants are 


SqA.-^ + a^A^ + agA^3 + a^X-^2 + a^A.-^ + a^ 

Ai = 

6A A-! 5 + 5BA. 1 1 ^ + ItCX-L 3 + + 2EA 1 


A 2 = 


a 0^2^ a i^2 a 3^2*~ a L^2 


6kXp + 52^2^ SLAg 2 + 2 EA 2 


Ao = 


SqA.^^ + a^A^ + + a 3^3^ + a lj^3 + a 5 


6k\p + 5BA. 3 ^ + iJCA 3 3 + 3LA 3 2 + 2E\ 3 



NACA TN 2129 


21 


^ *0^ ^ a^A>|^ O’ljAlj. ~t* 

^AXi,. 5 + 58X2^ + ilCX^ 3 + 3DX^ 2 + 2EX4 


A5. = — 
5 E 


A 6 


- l( at - a ^) 


The inverse Laplace transform of ^ is 


f = B ie XlSb + B 2 e X2Bb + B 3 e X36b + B^e^h + B 5^ + Bg 


( 8 ) 


where the .constants are 


Bi = 


hoX.5 + b-^X-j^ + hgX^ 3 + h^Xj 2 + hjjX^ + btj 
6AXL 5 + ^Xj 1 *- + IiCX-l 3 + 3BX-, 2 + 2EX^ 


BqX^^ t b 3_^2^" h^X^ 3 t b 3^2^ hjjXg ^ b^ 

Bg = ; 

t A Xp + 5B^2^ + l^Xg 3 + 3D*2 2 + 2E>2 


b 3 = 


^0^3"^ b i A 3 ^ **” hgX^ 3 + b^X^ 2 't b^X^ + b^ 
6AX 3 5 + 5BX 3 ^ + ^X 3 3 + 3DX 3 2 + 2EX 3 


BA = 


LqXa^ + h^X]^ + hgX^ 3 + h 3 X 4 2 + b^X^ + b^ 

6ax^5 + 5bxa^ + Ijcxa 3 + 3DXa 2 + sex^ 
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• B. 


5 = ^ 
p E 


B 6 - - b 5§) 


The inverse Laplace transform of p a is 


P = C^e + Cge^ 813 + Cj^e^ 813 + 


C 5 


vhere the constants are 


c qXi5 + c^X-^ ^ + c 2 ^ 1 ^ c^^i 2 t c j^X*^ 

6AX-l 5 + SBX^ + ilCX-L 3 + 3LX 1 2 + 2EX^ 


c 0^2^ c 1^2^ c 2^2~ ) C 3X 2 2 ^ J^X^ 

6AX2 5 + 5 BX 2 lt ' + itCXg 3 + 3 DX 2 2 + 2EX2 


CqX^ ; ~i~ CjXg^ ~f- CgXg 3 ~t~ CgX.j 2 c^X^ 

6ax 3 5 + 52 x 3 ^ + icx 3 3 + 3 dx 3 2 + 2EX 3 


n\- - C oH 5 + C 1^4^ + + + c ljAl). 

6 k\\p + 53 X 1 ^ + ICX4 3 + 3DX4 2 + 2EX4 


(9) 


c 5 = -= 


Cj}. 

E 
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The quantity Db0 can he obtained from equation (7) by differentiation 
as 


'■Db0 = Ai'e^ 1 ^ + A 2 ' e^ 8 ^ 3 + A 3 'e X 3 s b + A 4’e^ Sb + A^’ (10) 

where the constants are 

A ]_ 1 = XqAj. 

A 2 ' = X2A 2 

a 3’ = X 3 A 3 

A4' = X.4A4 

V - A 5 

The quantity hfcijr can be obtained from equation (8) by differentiation 
as 

= Bi'e^ 1 ^ + B 2 'e X2&b + B^e* 3815 + B^'e*^ + B 5 ’ (11) 

where the constants are 


B l' = X 1 B 1 


B ^ * — ^'2 B 2 


b 3 ' - x 3 b 3 

B^' = X4B4 


B^’ = B5 
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The quantity D^P can be obtained from equation (9) by differentiation 
as 


DbP = C 1 'a Xl ’ >b a a C 3 'e X 3 Sb + C| t 'a X, ‘ Bb 

where the constants are 


Ci' = X^ 


( 12 ) 


Cq 1 ~~ ^2p2 


c 3 ! 


C 1 — X4C4 


Collecting the equations of motion (equations (7) to (12)) for the case 
of four real roots gives 


i 

= 

+ A 2 e X2Sb a 

A 3 e X 3 8b a 

a/*’” 

+ Vb + 

y 

= B 1 e XlSb 

+ + 

Bja^ a 

Bi l e Xl ‘ Sb 

+ B 50b + 

p 

-•C 1 e X l s *> ‘ 

+ Cge^ a 

C 3 a X 3° b a 

C k e X ^ 

+ c 5 

i>b^ 

- A 1 '» XlSb 

a A 2 'a X2Sb 

aA 3 'e X 3 B » 

+ A^'e 

XX8b a A 5 1 

Dbt 

- B 1 ’e XlSb 

a B 2 'e X2Bb 

a B 3 , e > '3 B h 

+ B4V 

l 4®b + 

DbP 

= 0 1 'e XlB » 

+ C 2 'e X 2 Sb 

a C 3 'e X 3 Bb 

a cya X «> 




B 6 


\ (13) 


J 


Two Read. Roots and a Pair of Conjugate Complex Roots 

If a pair of conjugate conqplex roots X3 and X4 exists , the 
coefficients of the terms of 0, ty, and p (equation (13)) corresponding 
to X3 and X4 are conjugate complex. The complex number X 3 k can be 
written 

X 3 k = R k + Ifci 
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Thus, 

X 3 = Rj_ + I B i 

X^ — I?2 ^2^ 

X 3 3 = R 3 + I 3 1 
X 3 ^ = R^ + I 4 I 

X 3 5 = ^ + I 5 i 


The coefficients of the terms of 0 resulting from the complex roots 
are 


apRcj + ajR^ + a2R 3 + &3R2 + affix + + (apl^ + apIA + a2l3 + b.^2 + affix) 1 

A3 6 AR 5 + 5 BR 4 + 4CR 3 + 3DR 2 + 2ER-l + (6AI^ + 5 BI 4 + 4CI^ + 3DI 2 + 2EI 1 j i 


or after rationalizing 


A 3 - r a 


+ I A i 


and 


A b ~ r A " I A i 

Similarly, the coefficients of the terms of resulting from the 
complex roots are 


B3 = R B + I B i 

and 

B 4 = R B - Ip 1 

and the coefficients of the terms of (3 resulting from the complex roots 
are 

C 3 = Rq + I c i 


and 


c 4 = Be ” I C i 



2 6 
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The terms of <f) corresponding to the conjugate complex roots are 
A^e + A^e ^ &b = K A e Rl&b cosCl^s^, + m A ) 


where 


and 


K A = 2 \j B A 2 + J A 2 


cd a = tan - -* - — 

R a 

Similarly, these terms for ijr are 

B^e 3 + B^e*" 1 *^ 5 = Kpe Riab cos(l]_sb + eng) 

where 


and 


K B = 2 \ % 2 + 


i In 

cog = tan"- 1 - 

r B 


and for P are 

C^e ^ + C]^e ^ ^ = K^e^- 1 - 8 ^ cos(l^Sg + o^,) 


where 


and 





tan~-*- 


Rc 
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The final equations corresponding to two real and two conjugate complex 
roots are 


0 = A l e lSb + A^e^ 2 ^ + K A e RlSb cos (1-^ + ca^j + k^B h + Ag 


>l®b a. o X 2Sb ^ RlSb 


= B-je D + Bge d D + K B e 1 b 008^1-^ + + B^s^ + Bg 

0 = Cpe^ xSb + C2e^ 2Bb + K£.e Rl&b cos^I^&b + °^cj + C5 


Db0 - Ai'e^ + A 2 'e X 2 0 b + 1,^7^ ^ cos^i 


B b + “A + 


tan -1 — 1 + Ac 
Ri 




(1 


Db^ = + B 2 'e X2Sb + K^ 2 + I-^ e* 1 ** cos 


tan " 1 & 1 +B 5 


DbP = C^'e 1&b + C 2 'e ^ + Kq^|r ^_ 2 + I-j 2 e^ 18 * 3 cos^I-^sg + a^, + 


tan- 1 ^ 
Ri 


. y 
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Two Pairs of Conjugate Complex Roots 


If two pairs of conjugate complex roots, Xp, and X3, X4 
exist, another cosine term is introduced into equations (l4) in place 
of the exponentials so that 


<f> = K a ' e^- 1 - St> cos 
= K-g'e R l Sb cos 

p = Kq ' e Rb ^ cos 
Db0 = V^Rg' 2 + I 

Ka^I 2 + !l 2 


(il'Sb + “a') + K A eRlSb cos ( I l s b + “a) 


(ll ,s b + £ ^b') + K B e RlSb cosppsy, + tog) 


(il'Sb + (Pc') + Kce RlSb cos(lps b + a^j 


<2*1% coh^* 


Sg + + tan 


.1 V 

Ri ' 


e Rl0 b cosO^Bg + a> A + tan -1 ^ j + A 5 


+ A^sg + Ag 
+ B^Sg + Bg 


+ C^ 


>(15) 


= *b'\I*1 2 + *1 2 eRl % COB ^l'% + V + tan_1 5“T^ 


K- 


B 


\Jr^ + Ip 2 cos |lp&b + o>g + tan -1 


+ Be 


DbP = *c'\l R l' 2 + X l' 2 ^ % COB (*l'% + V + tan_1 


Kf 


;\Al 2 + Ip 2 e ElSb cos ^sg + + tan " 1 ^ 
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DISCUSSION 


The lateral motions of aircraft have been obtained by means of the 
Laplace transform. This analysis resulted in equations from which the 
free lateral motion of an aircraft can be calculated for any initial 
condition, or the forced motion can be calculated for any constant lateral 
force or moment applied at zero time. In general, the lateral forces and 
moments applied to the airplane by control movement or atmospheric tur- 
bulence are not constant but are arbitrary functions of time. After a 
solution has been obtained in terms of constant disturbing forces and 
moments, however, the solution for the arbitrary forces and moments can 
be obtained by Duhamel ' s integral as explained in references 4 and 5* 

Jhe nature of the motion indicated by equations (13) to (15) depends 
upon the form of the roots of the polynomial 

A = 0 

which is commonly referred to as the stability quartic. The roots of the 
quartic can take three forms - first, all four roots real; second, two 
real roots and a pair of conjugate complex roots; and third, two pairs of 
conjugate complex roots. In the case of the lateral motions of airplanes 
the first form almost never occurs; the second form is very common; and 
the third form occurs under rather rare conditions. The actual motions 
indicated by equations (13) to (15) can be seen to be composed, in general, 
of the sums of terms which are the amplitude coefficients (the A’s, B's, 
C's, and K's of equations (13) to (15)) modulated by exponential and cosine 
factors . 

All the classical stability concepts can be obtained from equations 
(13) to (15)- Because stability is concerned only with the free motion 
(motion due to Initial conditions) the forcing or disturbing quantities 
C 1 , C nc , and Cy c can be set equal to zero so that the amplitude coef- 
ficients A 5 , B 5 , and C 5 vanish. The variation of the amplitude of the 

motion with time, which determines the stability, is now dependent entirely 

X-i s-u XpBi. Xp^h ^k B b -^l 8 b 
upon the damping coefficients e , e , e , e , e , 

and e 1 ‘ 3 . The motion diminishes with time (stable) if Xp, X 2 , X 3 , 

X 4 ., Rp, and R^' are all negative. Thus, these criteria for stability 

are that all real roots of A = 0 and the real parts of all complex 
roots be negative. These criteria have been expressed in reference 7 in 
terms of the signs of the coefficients of the quartic A = 0 and the 
sign of Routh's discriminant which is written as 

R = BCD - AD 2 - KB 2 


(16) 
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These criteria in the present case can he expressed as follows: The 

necessary and sufficient conditions for the real roots and the real parts 
of the complex roots to he negative are that every coefficient of the 
quartic and R should he positive. 


If the motions contain oscillations , the periods of the oscillations 
in seconds are from equations (lL) and (15) 


P 


2rfb 

*1 V 


P' 


2jfc 

Il'V 


> (IT) 


and the times to damp to half-amplitude in seconds are 


, T l/2 = 


13 logeS 


R -jV 


T l/2 


h logg2 
R 1 , Y 



and the cycles to damp to half -amplitude are 


N l/2 


T l/2 = I 1 lQ ge 2 
P 2itRp 



Ij’ l°g e g 

2rtR l' 



( 18 ) 


( 19 ) 


APPLICATION 


The equations for the motion of an aircraft resulting from the 
analysis were used to calculate illustrative examples of certain free 
and forced motions of an experimental swept-wing airplane, a three-view 
drawing of which is shown in figure 2. The calculations were made hy 
use of the Bell Telephone Laboratories X-667Mf relay computer available 
at the Langley Laboratory. The calculations were based upon stability 
derivatives measured on a model of the experimental airplane in the 
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Langley stability tunnel and presented in reference 8. Motions were cal- 
culated for true airspeeds of 140 and 200 miles per hour under standard 
conditions. The stability derivatives, other related aerodynamic quan- 
tities, and the mass characteristics of the experimental airplane as used 
in the calculations are shown in table I. The coefficients of the sta- 
bility quartic (equation (4)) and value of Routh's discr iminant (equa- 
tion (l6)) were calculated as 



True 







C L 

airspeed 

(mph) 

A 

B 

C 

D 

E 

R 

0.693 

i4o 

26.19791 

10.18804 

3.021074 

0.6312249 

0.00235618 

5.8 

• 340 

200 

26.20030 

9.818377 

2.504971 

.4623735 

.00014875 

8.7 


The positive signs of all these quantities indicate complete stability 
of the lateral motion for both airspeeds. The coefficients of the quartic 
are such as to give two real roots and a pair of conjugate complex roots 
which are 



True 





airspeed 

(mph) 

X 1 

X 2 

• 

Rl ± llx 

0.693 

140 

-0.2802853 

-0.003603100 

-0.05249952 ± 0.28590791 

• 34o 

200 

-.2649690 

-.0003222716 

- . 05472583 ± 0 . 25197541 


Discussions of methods of obtaining the roots of the quartic can be found 
in references 9 to 12. 

The first or second powers of a multiplied by A in the denomi- 
nators of equations (3), (5)> and (6) introduce one or two zero roots, 
respectively, in addition to the roots given in the preceding table. 

These zero roots lead to the terms containing the amplitude coefficients 
A 5 , Ag, B 5 , Bg, and C 5 of equations (l4). For the experimental air- 
plane, the motion can be thought of as composed of three modes - the 
oscillatory mode resulting from the pair of conjugate complex roots, the 
rolling-subsidence mode resulting from the large negative real root 
and the spiral mode resulting from the small negative real root and the 
zero roots. Stability of the free spiral motion is indicated by the 
negative sign of the small real root Ag* 
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The period, time to damp to half -amplitude, and cycles to damp to 
half -amplitude of the oscillatory motion were calculated hy use of equa- 
tions (17) to (19) from the imaginary and real parts of the conjugate 
complex roots as 


C L 

True airspeed 
(mph) 

' P 
(sec) 

Tl/2 

(sec) 

N1/2 

0.693 

1*40 

3.60 

2.1 6 

0.60 

• 34o 

200 

2.86 

1.45 

■ .51 


The motions calculated fell into two categories which may he termed 
free and forced motions. 


Free Motions 

Free motions are those which exist following an initial condition 
and in the absence of any forcing function. The five possible initial 
conditions are f 0 , i}r 0 , p 0 , r Q , and p Q . Every free motion that the 

airplane is capable of executing can be obtained by superposition of the 
motions following these initial conditions taken separately. Figures 3 
to 6 show the calculated free motion following the initial conditions 0 q, 
Pq, tq, and pq for the experimental airplane in level flight according 
to equations (l*+). No airplane response to occurs when the angle of 

climb is zero. Table II gives the values of the amplitude coefficients 
(see equations (l*4-)) corresponding to the motions of figures 3 to 6. These 
figures show the total airplane motion and show the separate contributions 
to the motion by the rolling-subsidence and spiral modes when the motion 
resulting from these modes is appreciable. Figures 3 to 6 indicate that 
in the case of the experimental airplane the initial condition Pq pre- 
dominately excites the oscillatory mode of motion, tq excites both the 
oscillatory and spiral modes, and and po predominately excite the 

spiral mode. The rolling- subsidence mode appears for a very short period 
of time in the initial phases of any motion involving appreciable rolling 
velocity. The principal effects upon the motions of figures 3 "to 6 of 
increasing the airspeed from 1*40 to 200 miles per hour is to reduce the 
period as well as the time and cycles to damp to half- amplitude. 


Forced Motions 

Forced motions are those which exist during the action of forcing 
functions upon the airplane. Any forced motion of an airplane can be 
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built up by proper superposition (Duhaigel's integral, reference 4) of the 
motions following the application at zero time of constant values of the 
forcing functions C^ c , C nc , and CY c * Figures 7 and 8 show the cal- 
culated response according to equations (14) of the experimental airplane 
to the constant value 0.02 for. C| and C n applied at zero time. The 

response to a value of 0.02 for Cy c was also calculated but during a 
time period of 8 seconds was negligible compared with responses resulting 
from C 2 c and C n(; . For the experimental airplane, the value C^ = 0,02 

corresponds to a total aileron deflection of 21.0°, the value of C nc = 0,02 
corresponds to a rudder deflection of 13 . 7 °> and the value Cy c =0.02 
corresponds to a rudder deflection of 7*5°* The response to C is pre- 
dominately in the spiral mode of motion; whereas the response to C nc is 
predominately in the spiral and oscillatory modes. 

A large number of forced motions calculated for the experimental air- 
plane corresponding to various flight conditions are presented in refer- 
ence 8. These motions were built up by superposition of motions Buch as 
those of figures 7 and 8 following the application of the constant forcing 
functions Cj and C nc . A large number of comparisons are made in 

reference 8 between calculations and flight tests for a large variety of 
flight conditions. The agreement between calculated and flight motions 
is good and indicates the practicability of analyzing the dynamic lateral 
flying qualities of aircraft by use of the theory of lateral dynamics such 
as that herein developed if experimentally determined values of the aero- 
dynamic and mass parameters of the airplane are available. Reference 8 
also indicates rather insignificant effects upon the calculated motions 
that result from a consideration of slight nonlinearities which occur in 
certain aerodynamic parameters of the experimental airplane. 


CONCLUDING REMARKS 


The lateral motions of aircraft were determined by means of the 
Laplace transform which gave solutions expressed in terms of elementary 
functions for the free and forced motions. These equations permit the 
calculation of the free motion of an aircraft following any initial con- 
dition or the forced motion following the application of constant external 
forces and moments. These forced motions can be used to obtain the 
response to any arbitrary forcing function by means of Duhamel's integral. 
All the classical stability concepts can be deduced from these same 
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solutions largely "by inspection. These equations for the lateral motion 
were applied to the calculation of the lateral stability of a specific 
airplane and to the calculation of certain of its free and forced motions. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base^ Va., April 1950 
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APPENDIX 

EXPANSION OF THE COEFFICIENT DETERMINANTS 
The stability qoartic coefficients are 
A = (k X 2 K z 2 - %z 2 ) 

■B = 2^ 2 j2C Yp ^K xz 2 - K x 2 Kz 2 ) + K xz ^C 2r + C^ - ^X^Dj. “ ^Z^Ipj 

C = ^Pb 2 ( K X 2c np “ K XZ c Z p ) + ^X 2 ( C n r C Yp " C np°Y r ) + (^p^P 

• %%) + (%% ‘ C V*B + Cn p CY p ‘ Cn p%) + 

.. |^(%Cn r - C^) 

D = 2M b C L J^ 8 - 11 7 ( K XZ C 2p " + ( K XZ C np “ Kz^Z^jj + ^ZpPnp " 

CZpCnp^ + ^|pZ r CnpCYp + c Z p c np C T r + c ZpCn r CY p - c ZpC n pCy r - 

C Z p C %% - C Z r c np C Y p ) 

E - 7 (%S - %%) + % C *T - C *rS] 
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The coefficients appearing in the numerator of amplitude coeffi- 
cients for 0 are 

ao = ^o 8li b 3 ( K X 2K Z 2 “ k XZ 2 | 

a l = [2^ 2 Kxz(Cz r + C np ) + H^y^xz 2 - KX^Z 2 ) " 2 ^ 2 ( K X 2 °n r + 

Kz^Zp)] + Mo ^ 3 (% 2 K Z 2 " K XZ 2 )] 
a 2 = ^o[%Kx 2 (Cn r CY p - C^Cy^ + W^X^np - K XzCz p ) + tW^Zp^V 
C ip Cy P ) + ^XZ (CnpCYp - CnpCYp + CipCY r - Cz r CYp) + |w>(cZp c n r - 
c Z r c n p )| + (Hb0) o ^2CYp(KXZ 2 - ^X^Z 2 ) + 2hb 2 (KxzCz r - Kx^n*.) + 

(r>bt) 0 j2hb 2 (K Z 2 C Zr - KxzCnr) + P 0 [W^K Z 2 Czp - KxzCnp) + 

Vb 2 ( k Z 2c Zc “ K XZ C n c ) 
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a 3 = 0 ( 0 2^b tan 7C L (KxzCZp - Kj^C^ + Hb(c 2 pCn p - C 2 p Cnp) + 

( C l r %% + %% C Y r + % C *r% ~ %% Cy r " % Cl 3 r% " 

c z r c np c y p jj + Mo [^b 2 - K XZ c lp) + ^^(pn^Yp - 

C np C Y r ) + ^Xz( C lp°r r w C l r C Yp) + *0 j^Mb tan 7 C L ( K Z 2c 2p " 
K XZ C n p )] + (^)o[^Z 2 ( C Zp C Y r - C 2 r C Y p ) + V(%Z C n p “ %*%) + 

^XZ(V^ ' %%)] + P o[^( Ci rS" %%•)] + 

2P-b K Z S (% C Y c " %%) + 2|i b K Xz( C n c c Yp " C np C Y c ) + 

( C ^r Cn c " Cz c Cn r) 
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= 0o[§ tan 7C L ( c l p c np - %%)] + j^b tan ? C L (K E c l f 


Kx 2 ^)] + *o[§ r 


(Db*)o 


2Mt) tan 7Cjv, + 2^^C7 C T ,„ - C 7 „C.„ 1 + 


V Bp - U lp 


^c) 


^Ci c (CrL r CYp - CnpCY^ + | C^Ci^ - Ci r C Y p ) 

I <*cKS - %<=%)] 

*5 * rCl (%% - °I C S) 


The coefficients appearing in the numerator of amplitude coefficients 
for tJt are 


b o = J 8 ^ 3 (kA 2 “ K XZ 2 )] 

hi = t 0 [ 2 ^^Xz(cnp + C z r ) - ^Pb^p^X^Z 2 ~ Kxz 2 ) - 2 Pb 2 (Kz 2 Cz ] 

% 2 Cn r )] + (^b^of 8 ^ 3 ^ 2 ^ 2 " K XZ 2 )] 
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*2 - - K XZ%)] + %[^(VV % %) + 

MbK Z 2 fcl p Cl p - Oi p CY p ) + WjKxz (CnpC'fp - %% + r -lff?l T - Cl r CY p ) + 

- Kx2C Ip ^ + \ ^ ( c l p c ri r " c l r c n p )j + 

(l>b*)o[f^ 2 ( K XZ C n p - Kz 2 Cl p ) - i H*b 2 CY j! ^ x 2 K z 2 - Kxz 2 ]] + 
f*o[“"‘b 2 ( K X 2 c n [ , - K XZ G Zp)] + l *Mb 2 ^x 2 c n c * K XZ C Z C ) 
b 3 = ^o^C^KxSCng '- KxzClpJ + (Dbff) 0 [wx 2 (cnp(>f p - Cp^) + 

^XZ^Zp^p - %%)] + +o[^( c lp C n p - %%) - 2 Mb c L (*Z% - 
KxzS,) + 5( c Z r %% + 0 W*r + < W*p • 'Wr ‘ 

CZpC^CYp - Cz^CnpCxJ + (Db'Oo^Z^CljCYp - C I p CY p ) + 
t*b K Xz( C Hp C Y p - %%)] + ,! o[^b( C! p C Pp ‘ %%)] H 
W'( C V> - % c np) + ^XZ^l^Yp * C lp C X c ) + 

2 '% K X 2 ( c np C Y c - V^p) 
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H = ^O^L^lpCnp - Ci p Cnp| + (Db0) o ^PbCL^x 2 ^ - KxzCZpjJ + 
to^L^IpCnr " G Zr C ng| + (*>& ^RjCL^KxzCng - Kz^Zpjj + 
icz c (c np CY p - Cn p CXp) + |Cn c ( c 2 p c Yp ~ %%) + 

ic Yc (cz p C np - C ZpCnp ) 

= c l(cz c c iip - G Zp C n c ) 

The coefficients appearing in the numerator of amplitude coeffi- 
cients for 0 are 

c 0 = Po[8^ 3 (kx 2 k z 2 - kxz 2 J 

Cl = 0 o^ 2 c l (k x 2 Kz 2 - %z 2 )] + (Db0)o [s^ 2 Cy p (kx 2 Kz 2 - Kxz 2 )] + 

[W tan 7 Cl(kx%Z 2 - KXZ 2 )] + (^) 0 j 2 ^ 2 - *H*b) (k X%Z 2 - 

*4] + P 0 |2^ 2 (KXZCz r - Kx 2 Cn r ) + 2^ 2 ^CxzCn p - K Z 2 Cz p )J + 
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c 2 = ^o[i> c L^XZCz r - Kx^nj.) + *> cl ( ke % - **»,)] + 

(Bbf()o[|w : X 2 (cn p Cr r - Cn r CY p ) + SPb^x 2 (2K Z 2 C L - C n \ + 
2Wb%E ( C Ip - ^l) + iVxz(Cl r CY p - C lp C Ir jJ + 

+ 0 [ M b tan rc L ( K XZ C l r - C tl/-7 2 ) + Pb tan /C L - K 7 , 2 0 1{) )J + 

- S%) + ^(' lr % - %%) + 

4^ 2 tan 7 C l (k x % z 2 - K^ 2 ) + 21^2^2^ _ R^C^J + 

P 0^( C * p Cn r ' C2 r Cn p)] + H 2 (kxzCz c - Kx 2 ^) + ^z 2 (c Zc CY p - 

Cz p c Y c ) + ^ K X 2 (Cn c C Yr - V^ c ) + ^ K Xz(Cn p C Yc - C^CyJ + 

^XZ^Zr^c ~ c Z c c Yr ) 
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c 3 = ^o[| t! I.( C Jp C % - C J r %)] + Mjjyi- ta r ( K X 2 °n p - KXZC Ip ) + 

ZV>J K XZ c Z r - K^Cnj.y] + + 0 [| tan 7C L (c lp Cn r - Oj r Cn p ^j + 

(»b+)o[^ c L( K Z 2c l r - KjtzC^.) + 1^ tan 7C L (k^O^ - Kz 2 ^)] + 

^(cz p On c - c l c C Dp) + 2|i » CL ( Kz2Cl c ' K XZ c n c ) + 

21 % tan 7 ('l ( K x 2c n c - K XZ c l c ) + |( c ! r , C n p C Y r + c 2 p c n r °Y c + 

/ 1 

c Z r c ri c c Y p " c l r c n p c Y c - °Z p c n c c Y r " c Y c c n r c Y p ^ 

C4 = ^(Oj^ - Cj c C % ) + | tan 7C L (Cj^ - Oj^ 
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TABLE I 

AERODYNAMIC AND MASS CHARACTERISTICS OF AN 
EXPERIMENTAL SWEPT -WING AIRPLANE 


NACA TN 2129 


.b 

ll*o 

200 

w 

8700 

87OO 

s 

250 

250 

m 

270.2 

270.2 

P 

0.00238 

0.00238 

7 

0 

0 

b 

33.6 

33-6 

c L 

0.693 

O.3IK) 

a 

10.6 

4.8 

T) 

11.05 

5.25 

Pb 

13.51 

13.51 

V/b 

6.111 

8.730 

KX 2 

0.02329 

0.02219 

K z 2 

0.05932 

0.060 1*2 

K xz 

0.007316 

0.0035^ 

c h 

-0.0659 

- - . 

-0.0275 

c np 

0.100 

0.0975 

% 

-0.739 

-0.722 

% 

-0.325 

-0.31 

S 

-0.1 

-O.06 

% 

0.41 

0.3 

C 2 r 

0.12 

0.07 

C *r 

-0.280 

-0.280 

Cy r 

O.36 

0.1*0 


Miles -per -Hour units. 
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AMPLITUDE COEFFIdEHTB OF FHEE AUD FORCED MDTI033 OF EXPER3MEKEAL SHEPT-WIEO AIRPLANE 


(a) V => 140 piles per hoar. 


TnlfrlmT conditions 

end farcins Variables Bolllne-Baoaidence 
funationa node 


f 0.0*073926 

t Bi -.00222650 

fo P c x -.00131258 

P Ai' -.06978068 

r B,* .00381366 


Ap -.01808626 

B x .0097328* 
.00573756 
.30503*62 
-.01667089 


Ai -.2713271* 

B x ' .01*82880 

Ci .0087*177 
Aj.' .*6*73072 
-.02539885 


-.2127315* 
Bp .01162773 
Cx .00683*35 
Ax' .366*0206 
Bx' -.0199x603 


•333*235 
B X -.01931536 

Cx -.01138685 

Ax' —6053610* 
.03308*6* 


.07219731 

-.0039*581 

-.00232607 

-.12366306 

.00673858 


.00235150 

-.00012031 

-.00007576 

-.00*02776 

.00022013 


Ofl dilatory cede 



*A 

K B 

*C 


fXVBi 2 + n 2 

1 2 + z i 2 


K A 

Kc 

Ka\> 1 2 +i 1 2 
K B fl 2 + Il 2 



*A 

Kb 

Ko 

*a'Ai 2 +i i 2 

Krf/Hl 2 + Il 2 


Ka 

kb 

Kc 

K^l 2 + Ii 2 

KbK 2 + i i 2 


k a 

Kb 

Kc 

Ka^1 2+ . i 1 2 


0.05*0*332 A 2 

.0*0096*8 Bg 
.0*330260 C2 
.09600*16 Ag 1 
.07122*81 Bg 1 


.2*50096 a 2 
.1817706* Bg 
.19631*8* Cg 
.*332*085 Ag 
.3229020 Bg’ 


.02*12880 Ag 
.01790107 Bg 
.019333*0 Cg 
.0*266193 Ag* 
.03179911 Bg 1 


•33205361 Ag 
.2611839* Bg 
.28208*36 
.62538141 

.*6396533 


.07815380 Ag 
.05796158 Bg 
.06262090 Cg 
113883*29 Ag 1 
.1029990 S2 1 


.1935923 

.1*362*8 

.1551168 

.3*3902*0 

.25313879 


.003119*0 

.00231*25 

.002*99*3 

.0053*138 

. 00*11110 


0. *37*6*7 A5 
-3.038911 B3 
.01392006 
-.009632*9 
.066913*9 


-.02*58202 Aj 
.17076788 B5 
.00078222 C3 
.0005*129 A3 1 
—00376112 B5 1 


2817596* 

-1-9572863 

.00896556 

-.00620389 

.0*309623 


.*6209967 A 5 

-3-2100*13 B5 
.01*70395 
-01017*68 
.07067983 


A6 

0 

*6 

3.029296 



.025250*9 Aj 
-.175*060 Bj 
.000803*7 C5 
-00053599 A5 1 
.00386225 Bj 1 



*6 

-.0283600* 

*6 

.17*8510 


^NACA^d 
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TABLE II 

/KPLTTODE COEFFICIENTS 07 FREE ABD rOBCED >5DTIOH9 OF EXPERIKEBTAL SVEPT^WIBa AIHPIAHE - Co nclu d e d 


(b) V - 200 rrllf per boor. 


Till tl "1 ccnditiOTifl 
end forcing 
function* 

Variable* 

Rol 3 

node 

Ooclllatary code 



Spiral mode 




i 

A 1 

0.01197011 


0.01660600 

a 2 

0 . 6807 &X 

a 5 

B 

A6 

0 


+ 

Pi 

-.00027 171 

A B 

.02367936 

*2 

- 18.21875 ■ 

®5 

B ■ 

H 

1B.21629 

*> 

P 

Cl 

-.0003213 

Ac 

.02626556 

=2 

.0083529 

=5 


m 



P 

A l’ 

..02769022 

KX^Bl 2 + Il 2 

.03693086 

a 2* 

-•0013J277 

* 5 ' 

E 




r 

®1* 

.00069^6 

r B\ lR l 2+ Bl 2 

.05285335 

B2‘ 

.05126683 

* 5 ' 

Bl 

B 



i 

A 1 

-.10002626 

*A 

.13667276 

a 2 

-.00692530 

*5 

Bl 

m 



♦ 

Bl 

.00180871 

a b 

.19265072 

®2 

.18665768 

=5 


B6 

B 

Po 

P 

Cl 

.00270588 

Ac 

.19889500 

C2 

-.00000550 

C5 

E 

m 

E 


P 

A l* 

.23323073 

Dflsgsa 

.30270576 

a 2 > 

.00001386 

* 5 ’ 

B 




r 

Bl’ 

-.00618398 

*B(fSl 2 + Il 2 ' 

.63321690 

v 

-.00052525 

B 5 ‘ 


B 

■mu 


i 

Al 

-.2093808 

»* 

.00718425 

a 2 

.2112909 

*5 

Bl 

Afi 

0 


t 

% 

.00375619 

kb 

.01028182 

B2 

-8.006691 1 

=5 


B6 

7.995698 

Bo 

P 

P 

Cl 

Al* 

.00361936 

.10133861 

Ac 

A A \K 2 +'Il 2 

.01062610 

.01617213 

C2 

A2' 

.00367088 

-.00059656 

C 5 

AS' 

■ 

H 



r 

Bl’ 

-.00868879 

> J: 

.02316696 

B2' 

.02253050 

B 5 ' 

IB 




* 

A 1 

-.07670895 


.1518581 

a 2 

.1869736 

A5 


□ 

0 


* 

% 

.00137619 

A B 

.2173312 

®2 

- 7.085196 

B5 


»6 

7.097827 

r 0 

P 

Cl 

.00205853 

Ac 

.2266086 

C2 

.00326862 


1 

M 



P 

*1* 

.17766277 

Aa'/ r i 2 + i i 2 

• 36186086 

a 2* 

-.00052616 

*5 ■ 





r 

Bl' 

-.00318338 

Bi angina 

.68922668 

*2’ 

.01993765 

B 5 ' 


B 

BB9 


i 

Al 

.6567069 

Aa 

.03167098 

*2 

- 365-6037 

*5 

0 

A6 

365.1805 


* 

Bl 

-.00615719 

Ab 

.06503932 

»2 

13855.66 

”5 

6.657163 

B6 

-13855.50 

% 

p 

Cl 

-.01220331 

Ac 

.06656752 

=2 

- 6.351295 

P 5 

6.600000 

9 

— - 


P 

A l* 

-1.0518225 

.01886912 

Aa'| s 1 2+i 1 2 

.07086298 

A 2 I 

1.0286136 

* 5 ’ 

0 


— 


r 

B l’ 

Ab/Bi 2 + i 1 2 

.10138614 

*2' 


® 5 ’ 






0 


.03526760 

Aa 

.1276666 

a 2 

- 102.7051 

*5 

0 * 

*6 

102.5660 


♦ 


-.00063272 

Ab 

.1826785 

B2 

3892.267 

B5 

1.257163 

B6 

- 3892.093 


P 

P 

r 

l 

-.00096653 

-.08158060 

.00166360 

Ac 

A A fl 2 + x l 2 

Atf/Hl 2+ Il 2 

.1887956 

.28733535 

.61121969 

C2 

*2' 

B2' 

-1.786201 

.28895732 

- 10.950758 

1 

1.600000 

0 

10.976986 

1 



a 

A 1 

.00160829 

a a 

.00193011 

a 2 

. 05655717 . 

Hj 


A6 

-.05882356 

°*c 


B 

B 

-.00002526 

-.00003779 

-.00325766 

.00005866 

a b 

Ac 

A A ^1 2 + I 1 2 

a b ^i 2 + Ii 2 

.00276228 

.00285677 

.00636679 

.00621801 

B2 

C2 

*2' 

V 

- 2.163383 

.00098269 

-.00015915 

.00603139 

El 

1 

1 
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Figure 3.— 


ip-, radians 
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es per hour. 
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Figure 8.— 
at zero 
7=0°. 


Total motion 

Spiral mode 

Roilrg-sjpsbence mode 
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I' 8 

^.4 
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-.4 

I .8 
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n 
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■ 
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B 



fi 
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a 


g 

a 
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■ 
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n 
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■ 

p 

s 

a 

■ 

■ 

u 

■ 

p 

5= 

s 

5 


a 

■ 

■ 

■ 

u 

u 

n 

a 
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■ 

■ 

■ 

■ 

■ 

■ 

■ 

u 
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2 3 4 5 6 0 

Time, seconds 
(a)V»m3es per hour. 


2 3 4 5 6 

Tme, seconds 
(b)V*200 mBes per hour. 


Response of experimental ' swept-wing airplane to the applicatior 
time of a constant yawing-moment coefficient. C ric = 0.02; 
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